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Abstract 

We consider the fragmentation process with mass loss and discuss self-similar properties of the arising 
structure both in time and space, focusing on dimensional analysis. This exhibits a spectrum of mass 
exponents 6, whose exact numerical values are given for which or t^' has the dimension of particle 
size distribution function tp{x,t), where z is the kinetic exponent. We obtained conditions for which the 
scaling and fragmentation process altogether break down and give explicit scaling solution for special 
case. Finally, we identify a new class of fractals ranging from random to non-random and show that the 
fractal dimension increases with increasing order and a transition to strictly self-similar pattern occurs 
when randomness completely ceases. 

PACS number(s): 05.45.Df, 46.50. -fa, 62.20.Mk 

The kinetics of the irreversible and sequential breakup of particles occurs in a variety of physical processes 
and has important applications in science and technology. These include erosion iQ , grinding and crushing of 
soHdsQ, polymer degradation and fiber length reduction jsj, breakup of liquid droplets [Q, etc. to name just 
a few. In recent years, there has been an increasing interest in studying fragmentation, allowing variations 
that increase the flexibility of the theory in matching such conditions of real phenomena such as extension to 
higher dimension |^ , agglomerate erosion Q , mass loss Q , volume change and fragmentation-annihilation 
The kinetic equation approach of fragmentation is described by a linear integro-differential equation 
where mass or size is the only dynamical quantity. This theoretical approach is mean field in character since 
fluctuations are ignored altogether. 

In this article, we study the kinetics of fragmentation with continuous mass loss, an interesting variant of 
the classical fragmentation equation introduced by Edwards et al. ^J. This is relevant in all fragmentation 
processes mentioned earlier where mass loss might occur due to evaporation, oxidation, sublimation, dissolu- 
tion, melting etc., or in Yule- Furry process of cosmic shower theory where energy loss occurs due to collision 
. Although there exist a series of paper devoted to this problem, the explicit scaling solution with exact 
numerical exponents, the geometric properties of the arising pattern and the interplay between fragmentation 
and mass loss remain unexplored. In addition to obtaining these, we give an alternative interpretation of the 
existing scaling theory by considering dimensional analysis, and we give an explicit bound to the exponent 
of the breakup rate for which the kinetic equation fails. A strong motivation for the present investigation 
came from the desire to know how fractal geometry of the resulting object changes with the degree of order, 
which is quantified by the global exponent and is typically known as fractal dimension. 
The evolution of the particle size distribution function -ip^x^t) for fragmentation with mass loss in one 
dimension is 

*^ = t) f F{y, X - y)dy + 2 / dy^{y, t)F{x, y ~ x) + ■^{m{x)^{x, t)), (1) 
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where F(x, y) is the breakup kernel describing the rules and the rate at which a particle of size (x + y) breaks 
into sizes x and y. Equation (1) describes a process whereby cuts are equivalent to seeds being sown on 
the fragmenting objects, thus producing two new segments. This immediately creates two more new ends 
belonging to the two different, newly created fragments; in doing so, fragments start losing their masses 



immediately from either end until they encounter another seed or become dustlike, thereby stopping the 
loss of their masses. Equation (1) also describes the sequential deposition of point-size particles that grow 
once deposited successfully and stop growing upon collision with another point or growing particle. In other 
words, the present model describes nucleation and growth when 'ip{x,t) describes the gap size distribution of 
size X at time t or how space is covered by growing objects. This occurs in a number of natural phenomena, 
including phase separation, wetting, droplet growth, and growth of breath figures. Recently, a variant of the 
present model was considered in pO| in which the deposition of growing rods instead of growing seeds was 
addressed. 

We consider the breakup kernel to be F{x,y) = {xy)^{x + yY^^ ^ for which the breakup rate a{x) = 
F{y, X — y)dy — px^^'^''^\ where p = ^r|f^_j^2) ' ^^"^ homogeneity index A(/3, a) — 2(3 + a, and > —1. 
To comply with the present choice of breakup rate a(x), it is essential to consider a similar power law form 
for m{x), hence we choose m{x) ~ mx'^ , with m a positive real constant. However, a{x) being the quantity 
describing the rate at which particles are fragmenting, x^^^^'"^"^ = t{x) must bear the dimension of time, and 
this put a strong constraint on the exponent 7. That is, the dimensional consistency requires 7 — A(/3, cr) + 1. 
This dimensional consistency has been ignored in all previous studies instead 7 < A(/3, c) + 1 was 
identified as the recession regime and 7 > A(/3, cr) + 1 as the fragmentation regime. A closer look at the 
rate equation reveals that x and t are inextricably intertwined via the dimensional consistency, and hence 
the system becomes stochastic in nature. So, it is obvious that either of the two can be taken to be an 
independent parameter when the other one is expressible in terms of this. For example, if x is chosen to be 
the independent parameter, then using the fact that the dimension of a physical quantity is always expressed 
as power monomial, we can define the dimensionlcss quantity ^ = It is obvious that the dimension of 
the governed quantity ^'{x, t) as well can be expressed in terms of x alone, and hence we can define another 
dimensionlcss quantity 11 = ^^jf' ~ x^iIj{x, x"^) = F{x, ^). Since ^ and 11 are dimensionlcss quantities upon 
transition from one system of units of measurement to another inside a given class, their numerical values 
must remain unchanged meaning ^ = 0. This implies that 11 is independent of x and can be completely 
expressed in terms of ^ alone. Thus we can define 11 — to enable us to write the spatial scaling ansatz 
ip{x,t) ~ x~^^{(^). Had we chosen time t to be the independent quantity, a similar argument which would 
lead us to write the temporal scaling ansatz ^{x^t) ^ t^^(l){ri) with 77 = p^. We know that has the 

dimension of time, therefore t ^(f ") = S{t) must have the dimension of x. This gives us 2; = jj^-^, which 
is known as the kinetic exponent since 6{t) describes the mean or typical cluster size and a = — A(/3, cr). 
Inserting the temporal scaling ansatz into Eq. (1) and requiring scaling to exist would give the same result 
for kinetic exponent z. However, dimensional consideration proved to be very instructive and time saving, 
yet rich in physics. 

Note that the exponent 9 takes the value for which x~^ and t^^ have a dimension of ?/'(x,t), and hence it is 
called the mass exponent. The existence of scaling or a self-similar solution actually means that we can choose 
self-similar coordinates (or -^^w) and x/5{t) (or t/T{x)) such that their plots for any initial condition 
collapse into one single curve. It is very instructive to note that at A(/3, cr) = 0, the governing parameters 
t and the governed parameter -0 all lose their dimensional characters. As a result, the system loses its 
stochastic nature. This means that one can no longer define self-similar coordinates at A(/3, cr) = 0, which is 
conceptually very important and crucial for scaling to exist, and hence scaling at A(/?, cr) < breaks down. 
Note the inherent properties of the fragmentation process, which are the typical or mean cluster size 5{t), 
must decrease as time proceeds, and the number density must be an increasing function of time. In order to 
be so, the system must be governed by some conservation laws. However, this is not true as the homogeneity 
index A(/3, cr) < 0, in which case 5{t) becomes an increasing function of time. This simply goes against the 
principle of kinetics of fragmentation. Combining all this, we argue that the system fails to show scaling 
not only at A(/3, cr) < 0, as was reported by Cheng and Redner but also at A(/?, cr) = 0. Many authors 
noticed further anomalous behaviour in this regime, e.g., lack of conservation of mass and an absence of 
self-averaging |]l3j. Due to this anomalous behaviour, this regime was termed a transition to shattering jl^ . 
These heurestic arguments actually imply a possible failure to describe physically meaningful fragmentation 
process by the breakup kernel, for which the homogeneity index A(/3, cr) < and the shattering is actually 
an articulate term for this regime. The above discussion is true for Eq. (1) even when the mass loss term is 
absent. Therefore, we shall restrict the rest of the discussion to the regime A(/3, a) > with 3?(/9) > —1. 
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We now turn to find the mass exponent 0, which can only be found if the system follows some conservation 
laws. For example, for pure fragmentation (m = 0), the mass or size of the system is a conserved quantity 
and gives 9 = 2. Defining the n*'* moment Mn{t) = x""ip{x, t)dx and combining it with the rate equation 
(1) for the present choice of F{x, y) and ■m{x) yields 

-ItT ^ -[r(2/3 + 2) r(n + 2/3 + 2) + ^n]M^^,,^^it). (2) 

Note that the number density Mo(i) evolves in the same fashion as it would in the absence of the mass 
loss term. This means that particles keep losing their mass in a continuous manner by some mechanisms 
that do not alter the number density. The interesting feature of the above equation is that for m > 0, 
there are infinitely many n = Df(f3,m) values for which A/£)^ (^ „i)(i)s are conserved quantities. However, 
for TO = 0, there is only one conserved quantity Mi(t), i.e., size or mass of the system, and this does not 
depend on /S. The dynamics of the system is governed by conservation laws and, as a consequence, the 
system shows scale invariance. These conserved quantities in fact are the intrinsic agent responsible for 
tuning the numerical value of the mass exponent 9, and introduce universality to the process. However, 
at A < there is no evidence that there exists any conservation law for which the mean cluster size can 
be a decreasing function of time while the mean number density can behave in the opposite way. We can 
find the numerical value of Df{(3,m) by searching for the positive and real root of the equation obtained 
by setting the term in the bracket of the Eq. (2) equal to zero, which is polynomial in n of a degree 
determined by the /3 value. Substituting the temporal scaling ansatz into the definition of Mn(t) gives 
Mn{t) ~ 77"(/)(7y)d77, and demanding MDf{(3,m) be a conserved quantity immediately gives 

9 = {1 + Df{l3, to)), which clearly depends on /3 and to if and only if to > 0. Owing to the random nature 
of the process and due to the presence of mass loss term, it is clear that when the process continues ad 
infinitum, it creates a distribution of points (dust) along a line at an extreme late stage that is different 
from any known set To measure the size of the set created in the long-time limit, we define a line 

segment S{t) = — which is typical cluster size. We can count the number of such segments 

needed to cover the set, and in the limit S — > (i.e., t — > oo), the number N{S) will simply measure the 
set and appear to scale as N{S) ~ ^-Dsifi^"^) , The exponent Df{[3, to) is known as the Hausdorff-Besicovitch 
Ip^ dimension or the fractal dimension of the arising pattern. 

To get a physical picture of the role played by to, we set /? = for the time being, for which the polynomial 
equation becomes quadratic in n and the real positive root is Df{m) = —^{l + l/m) + ^y/[l + I/to)^ + 4/to 
when the second root is D = —{Df{m) + 1 + 1/to). Therefore, the exponent 9 is also function of to. The 
expression for Df{m) reveals that as value of to increases, the fractal dimension decreases very sharply and 
in the limit to — > oo, Df{m) — > 0. This means that as to increases, the size of the corresponding arising 
set decreases sharply due to fast disappearance of its member, whereas, as to — > 0, Df{m) — > 1, that is, 
we recover the full set (pure fragmentation) that describes a line. On the other hand, had we kept to fixed 
and let p decrease, the effect would have been the same as we observed for increasing to with p = I (i.e. 
(3 = 0). Thus, it is the ratio between to and p that matters rather rather than their individual increases 
or decreases. To give further physical picture of what these results mean, we define mass length relation 
for the object as Mq ^ and Me ~ S"^ for the space where the object is being embedded, where d 

describes the Euclidean space. The density of the property of the object p then scales as p ^ ^Df{m)-d^ 
is thus clear that for a given class of set created by a specific rule, when Df{m) decreases it means that it 
is increasingly moving away from d and hence more and more members are removed from the full set. This 
in turn creates increasingly ramified or stringy objects, since Df{m) = d describes the compact object with 
uniform density. Therefore, this shows that increasing m/p ratio means that the mass loss process becomes 
stronger than the fragmentation process and vice versa. 

We now attempt to find the spatial scaling solution for $(C). Note that the dimension of the arising pattern 
is independent of cr and consequently independent of how fast or slow the system performs the process. 
Therefore, we can set a = \ without risk of missing any physics, but it certainly simplifies our calculation. 
Substituting the spatial scaling ansatz into the rate equation (1) for F(x, y) = \ and m{x) = mx^ and 
differentiating it with respect to ^, transforms the partial integro-differential equation into an ordinary 
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differential equation, 

^(1 - mO^"{0 + {(1 C[2m(2 -9)- 1]}<I>'(0 - [m(2 - 9){1 - 6) - {3 - 0)]<I>(e) = 0. (3) 

For m — 1, this is hypergeometric differential equation [|l7| whose only physically acceptable linearly indepen- 
dent solutions are 2^1(1, -(l-h2D/);-D/;0 and ^^1+ " 2i^i(2-hD/, --D/; 2-hD/; where D/ = 0.414213. 
From these exact solutions for spatial scaling function, we can obtain the asymptotic temporal scaling func- 
tion e~^i^ which satisfies the condition — > as ^ — > oo. 

We now attempt to see the role of f3 on the system. To judge its role, it is clear from the previous discussion 
that we ought to give equal weight to all the terms in Eq. (1) so that each of them can compete on an equal 
footing. This can be done if we set m — p — ^(2/3+2) that the relative strength between fragmentation 
and the mass loss process stays the same as the value of P increases. This is a very crucial point to be 
emphasized. We can obtain the fractal dimension for different values of which is simply the real positive 
root of the polynomial equation in n of degree (3. A detailed survey reveals that the fractal dimension 
increases monotonically with increasing (3. To find the fractal dimension in the limit /3 — > 00, we can use 
the Stirling's approximation in the polynomial equation to obtain ln[n + l] = (l — n) ln[2] when n = 0.4569997 
solves this equation. In order to give a physical picture of the role of f3 in the limit /3 — > 00, we consider 
the following model: F{x, y) = {x + y)^5{x — y). This model describes that cuts are only allowed to be in 
the middle in order to produce two fragments of equal size at each time event. This makes a{x) = , so 
we need to choose m{x) — ^a;''"*"^, {m — ^ gives the same weight as for the fragmentation process). Then 
the rate equation for Af„(t) becomes M^{t) = —[{n+ l)/2 — 2~"']Mn+-){t). As before, we set the numerical 
factor on the right-hand side of this equation equal to zero and then take the natural logarithm on both sides 
to obtain the n value for which M„(t) is time independent. In doing so, we arrive at the same functional 
equation for n as we found for (3 — > 00. This shows that the kernel F{x,y) — {xy)^{x + yY~^ behaves 
exactly in the same fashion as for F{x, y) — {x + y)^6{x — y). We thus find that in the limit f3 — > 00, the 
resulting distribution of points is a set with fractal dimension Df = 0.4569997, which is a strictly self-similar 
fractal as randomness is seized by dividing fragments into equal pieces. We are now in a position to give a 
physical picture of the role played by (3. First of all, the process with /3 = that describes the frequency 
curve of placing cuts about the size of the fragmenting particles is Poissonian in nature. Consequently, the 
system enjoys the maximum randomness and the corresponding fractal dimension is Df = 0.414213. For 
/3 > 0, the frequency curve of placing cuts about the size of the fragmenting particles is Gaussian in nature, 
meaning that as the value of f3 increases, the particles are increasingly more likely to break in the middle than 
on either end. That is, as f3 increases, the variance decreases in such a manner that in the limit f3 — > cx), 
the variance of the frequency curve becomes infinitely narrow, meaning a J-function distribution for which 
the fragments are broken into two equal pieces. Therefore, there is a spectrum of fractal dimensions between 
/? — > when Df = 0.414213 and f3 — > 00 when Df = 0.4569997. A detailed numerical survey that we do 
not present here, confirms that the fractal dimension increases monotonically with (3 and reaches a constant 
value when (3 — > 00. The previous density-dimension relation implies that increasing P vis-a-vis increasing 
order also means that the system loses less and less mass, and this happens despite the fact that now the 
^ ratio stays the same. This shows that there exists an interplay between fragmentation and the mass loss 
process that can be tuned either by changing the ratio of m and p, which is obvious of course, or by changing 
the degree of order alone, which is indeed a nontrivial result. 

In summary, we have identified a new set with a wide range of subsets produced by tuning the degree of 
randomness only. The process starts with an initiator of unit interval [0.1], and the generator divides the 
interval into two pieces, and deleting some parts from either side of both pieces at each time step. The 
amount of parts to be deleted is determined by the parameter that controls the intensity of randomness. 
We quantified the size of the resulting set obtained in this way by fractal dimension and showed that the 
fractal dimension increases with increasing order and reaches its maximum value when the pattern described 
by the set is perfectly ordered, which is contrary to some recently found results |l^. We also discussed the 
scaling theory of the process, emphasizing dimensional analysis, and we showed that the shattering is in fact 
is an articulate term whereby the equation fails to describe physically meaningful fragmentation process. 
We have also shown that the interplay between fragmentation and mass-loss arises not only from the ratio 
of their strengths determined by their respective numerical coefficients, but also by the degree of order. 
We give exact numerical value of the mass exponent, which has never been reported, and we obtained the 
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explicit scaling function for special case of interest. Finally, we argue on the basis of our findings that fractal 
dimension, degree of order and the extent of raniifications of the arising pattern are interconnected. 
The author is grateful to R. M. Ziff for sending valuable comments and acknowledges inspiring correspon- 
dence with P. L. Krapivsky. MKH acknowledges the Alexander von Humboldt Foundation for granting the 
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